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In Landau levels N > 1, the ground state of the two-dimensional electron gas (2DEG) in a 
perpendicular magnetic field evolves from a Wigner crystal for small filling v* of the partially filled 
Landau level, into a succession of bubble states with increasing number of guiding centers per 
bubble as u* increases, to a modulated stripe state near u* = 0.5. In this work, we show that these 
first-order phase transitions between the bubble states lead to measurable discontinuities in several 
physical quantities such as the density of states and the magnetization of the 2DEG. We discuss 
in detail the behavior of the collective excitations of the bubble states and show that their spectra 
have higher-energy modes besides the pinned phonon mode. The frequencies of these modes, at 
small wavevector k, have a discontinuous evolution as a function of filling factor that should be 
measurable in, for example, microwave absorption experiments. 

PACS numbers: 73.43.-f, 73.20.Qt, 73.21.-b 

I. INTRODUCTION 

It is now well established both theoreticallyi and experimentally 2 , 3 that the two-dimensional electron gas (2DEG) 
ground state near half filling in the higher Landau levels (N > 1) is the quantum Hall stripe state. Transport 
experiments have shown that this stripe state has an highly anisotropic longitudinal conductivity i.e. a yy » a xx 
(where x is the direction perpendicular to the stripes) and a Hall conductivity a xy that is not quantized. Away from 
half filling, more precisely at partial filling factors around i/* « 1/4 and v* w 3/4 in Landau level N = 2, the isotropy 
of the longitudinal conductivity is restored and a minima appears in the diagonal resistance. Concomitantly, the Hall 
conductivity is quantized, but at a value equal to that of the nearest integer quantum Hall effect plateaus. These 
other ground states around v* w 1/4 and v* w 3/4 in N = 2 have been called reentrant integer quantum hall states 
(RIQHS) 4 . Studies by Koulakov, Fogler, and Shklovskii and Moessner and Chalker 1 suggest that the RIQHS are due 
to the formation of bubble states. These states can be described as triangular Wigner crystals of electron clusters. 
Each cluster (or "bubble") contains a fixed number M of electrons such that, within them, the local filling factor is 
equal to one. The bubble states, being pinned by disorder, should be insulating in contrast to the stripe states. A 
review of these novel ground states in higher Landau level is given in Ref. H 

The ground state of the 2DEG changes discontinuously as the partial filling factor v* is changed. In the Hartree- 
Fock approximation (HFA), the ground state is a Wigner crystal at small filling v* or, equivalently, a bubble state 
with one guiding center per bubble. As the partial filling is increased, bubble states with increasing number M of 
electrons per bubble are stabilized so that the 2DEG evolves from the Wigner crystal, through a succession of bubble 
states, into finally the stripe state near half filling. In Landau level N, it can be shown that (in the HFA) the last 
bubble state has M = N+l. The optimal number of guiding centers in a bubble is well approximated by the formula 
M = 3v*N corresponding to an average separation between the bubbles equal to a = 3.3i? c where R c = ^2N + It is 
the cyclotron radius and t = ^Hc/eB is the magnetic length. This sequence of phase transitions in higher Landau 
levels has been verified by a number of authors. By comparing the energy of the Laughlin liquid at v* = 1/3 and 
v* = 1/5 with that of the bubble states in several Landau levels, Fogler and Koulakov^. showed, within the HFA, that 
the bubble states are lower in energy than the Laughlin liquid at filling factor v* — 1/3 for N > 2 and at filling factor 
v* = 1/5 for N > 3. Rezayi et ali used exact diagonalization studies on finite size systems to show the absence 
of fractional quantum hall states (FHQS) in higher Landau levels and a tendancy to charge density wave (CDW) 
ordering. Recent calculations by Yoshioka and Shibata— using the density matrix renormalization group method 
(DMRG) basically confirms the Hartree-Fock scenario with some minor corrections. For N — 2 and N = 3, the 
DMRG gives a wider region of stability for the stripe phase than in the HFA. As a result, the last bubble state (that 
with M = N + 1) is absent from the phase diagram. The transitions between different bubble states occur at slightly 
different filling factors in the two approaches. Both calculations, however, give first order phase transitions between 
the bubble states and between the bubble and stripe states. The possibility of a continuous phase transition involving 
the deformation of the bubbles from a isotropic shape to a more elliptical shape as the stripe state is approached was 
investigated (with a negative result) by Ren et al£ within the HFA. These authors also claimed^ that the v* = 1/6 
composite fermion state has lower energy than the corresponding bubble state in N = 2 and that this liquid state 
could be an intermediate state between the Wigner crystal and the M = 2 bubble state. (See, however, Ref. 0.) 
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In any case, the difference in the cohesive energy of the various possible ground states at any given filling factor is 
usually quite small and the possibility that quantum fluctuations, disorder, finite extension of the wave function in 
the third dimension, screening corrections, etc. modify the above scenario cannot be ruled out. 

The existence of the bubble state near filling factors v* sw 1/4 and v* 3/4 in N = 2, was first evidenced 
by transport measurements showing its insulating behavioi— and its nonlinear I — V properties^. In more recent 
experiments, Lewis and ali^i and Ye et alH measured the diagonal conductivity o xx of the 2DEG in the microwave 
regime and found sharp resonances in -ft [cr xx ] at some frequency. These resonances are strongest near c' « 1/4 and 
v* w 3/4 in level N = 2. They appear for a range of partial filling factors around i/*«l/4 and v* sa 3/4 and only for 
temperatures below 0.1 K. It is believed that these resonances are due to the excitation of the pinning mode of the 
bubble states realized at these filling factors. The experiments show no hint of discontinuous transitions among the 
different bubble states. It is possible that thermal or disorder effects may smear out any sharp effects on the pinning 
mode. 

In light of these microwave experiments, it is natural to ask how one could observe a transition between bubble 
states, and also how one could distinguish between bubble states of different M numbers. In this paper, we study how 
the behavior of several measurable physical properties of the bubble states evolve as the filling factor, i/* , is increased 
from the Wigner crystal state at low filling to the stripe state at half filling. We show that the series of transitions 
between bubble states lead to oscillations and discontinuity in the orbital magnetization and susceptibility and to 
features in the one-particle density of states that are directly related to the number of electrons per bubble. Working 
in the time-dependent Hartree-Fock approximation (TDHFA), we also compute the dispersion relation of the phonon 
as well as those of other higher-energy modes that correspond to excitations localized on the bubbles. These localized 
excitations are interesting because (in contrast to the phonon mode) they are gapped and probably less sensitive to 
the disordering potential. Their dispersion can be computed, in a first approximation, by neglecting disorder which 
is difficult to take into account in these charge density wave states. We thus apply the TDHFA to show that, within 
the range of frequencies accessible in microwave experiments, the first order transitions between the bubble states can 
be tracked by looking at the discontinuous behavior of the first few high-energy modes as a function of partial filling 
factor. 

The outline of the paper is as follows. In Sec. II, we briefly describe the Hartree-Fock procedure to compute the 
order parameters of the charge density wave states as well as the orbital magnetization and the one-particle density 
of states. We then introduce the time-dependent Hartree-Fock approximation and describe how we get the dispersion 
relation of the collective excitations in these ordered states. In Sec. Ill, we discuss our numerical results for these 
quantities. We conclude in Sec. IV, by discussing how our results could be related to the experiments of Lewis et al. 
and Ye et al.. 



II. HARTREE-FOCK FORMALISM 



A. Order parameters 

The bubble phase is a crystal state that can be described by the set of average values {(n (G))} where (n (G)) is 
the ground-state average of the density operator and G is a reciprocal lattice vector of the crystal. It is convenient, 
in the Landau gauge, to define a related density operator by the expression 

P(G) = -Lj2 e - iG ° X+iG * Gyll/2 4,x%,z- Gsii) (1) 
(p x 

where c N x creates an electron in Landau level N with guiding center quantum number X and N v is the Landau 
level degeneracy. This new density operator is related to the real density by the expression 

n(G) = N v F NtN {G)p(G), (2) 

where Fn,n{G) is a form factor for the N th Landau level which is given by 

F N , N {G) = e- G ^"L%{^fy (3) 

and LP N (x) is a generalized Laguerre polynomial. At a semiclassical level, the average (p(G)) can be viewed as a 
Fourier transform of a "guiding center density" of cyclotron orbits. 

In this paper, we make the usual approximations of neglecting any Landau level mixing and consider that the filled 
Landau levels are inert. It follows that the filled levels do not enter in our calculation and we will, whenever possible, 
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drop the level index N from now on to simplify the notation. Moreover, we assume that the partially filled level is 
completely spin polarized. 

The average values (p(G)) are obtained by computing the single-particle Green's function 

G(X,X',t) = -(Tc x (t)c* x ,(P)), (4) 

whose Fourier transform we define as 

G ( G ' T ) = W E e-^( x+x ')s x , x ,_ Gyll G(X,X',r), (5) 

* X,X' 

so that 

(p(G))=G(G,r = 0-). (6) 

We compute the single-particle Green's function of Eq. (J5J by solving numerically the Hartree-Fock equation of 
motion 

[iw n + /i] G (G,w n ) - Y, U ( G - G ') e* GxG ' £2/2 G (G'.iwn) = £ G ,o, (7) 

G' 

where u) n is a fcrmionic Matsubara frequency and the Hartree-Fock potential U (q) is given by 

U (q) = [H NtN (q) - Xjv^ (q)] (p (-q)) • (8) 
In Eq. (JSJ), the Hartree and Fock interactions are given by 

fa 2 f 



e 2 \ 1 



N 
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X^(q) = ^v^jr°°<i!ce- Ba [ZiS r (x 2 )] a Jo(v^r^), 



(9) 
(10) 



where re is the dielectric constant of the host semiconductor and Jo (x) is the Bessel function of order zero. 
The electronic density in real space is obtained by the relation 



< n W) = 2^2 E e+2G ' r ^ (G) <p (G)) . (11) 

G 

B. Hartree-Fock energy 

We solve the Hartree-Fock equation of motion (Eq. Q) numerically by using an iterative approach that was 
described in detail in Refs. • Once the order parameters{(p (G))} are found, the Hartree-Fock energy per particle 
in the partially filled level can be written as 

if 1 = 2^ E l H *» ( G ) C 1 - *6.o) - Xn^n (G)] I (p (G)> | 2 , (12) 

where i/* is the filling factor of Landau level N and N s is the number of electrons in the partially filled level. The 
total filling factor of the two-dimensional electron gas is v = 2N + v* . 

C. Density of states 

In our iterative approach, the single-particle Green's function is computed by first finding the eigenvalues and 
eigenvectors of the interaction matrix 

F G , G , = U(G-G')e lGxG ' 12 / 2 . (13) 
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This matrix is hermitian and can be diagonalized with the unitary transformation 

F = VDV\ 

where V is the matrix of the eigenvectors of F and Dij = djSij is the diagonal matrix of the eigenvalues of F. The 
Green's function is then readily given by 



G(G,W n ) = V: 
* J 1 



iw„ + (p — dj) /H 



(14) 



The density of states is defined as 



9{v 



) = ~J dx$S[G R (r,r, 



U3) 



(15) 



where G R is the retarded single-particle Green's function which can be computed from the eigenvalues and eigenvectors 
found above. We then have 



g(uj) = ^ 

IT 



\Vg=oj\ 



to + iS — dj/H 



(16) 



D. Collective modes 



As shown in Refs. 0, the order parameters {(p(G))} can also be used to compute the density-density response 
function xg.g' (k, r) defined by 



(k, r) - -N v (Tp(k + G, r) p (-k - G', 0)) , 



(17) 



where p = p — (p). The collective modes of the bubble state are found from the poles of this response function. By 
following the poles with non-vanishing weight as the wavevector k is varied in the Brillouin zone of the reciprocal 
lattice, we get the dispersion relation of the phonon and higher-energy collective modes. 

The equation of motion for the density response in the TDHFA is given by the matrix equation 



\}^nS G ,G' - M G ,G» (k)] XG",g> (k, iQ n ) = Bg,g< (k) , 



G" 



where Q n is a Matsubara bosonic frequency and the matrices Mg.g arid Bgc are defined by 



M g ,g' (k) 



■ 2i (^)<"( G - G ')> 



x sin 



(G x G') 



[H N . N (G - G') - X N . N (G - G') - H N , N (G') + X N , N (G')] , 



(18) 



(19) 



and 



Bg.c (k) = 2«sin 



(G x G')l 2 



(p(G-G')) 



respectively. 

To solve for XG,G' (k, iCl n ), we diagonalize the matrix Mq,g" (k) by the transformation 



(20) 



M = WTW~\ 



(21) 



where W is the matrix of the eigenvectors of M and Tjj = tjSij is the diagonal matrix of its eigenvalues. The analytic 
continuation of XG,G' (k, iQ n ) is given by 



Xg,c (k,a;) = ■ 



w Gd {k) \wiky 1 



3>k 



Bk.G' (k) 



lu + iS — tj (k) 



(22) 
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The jth eigenvector in the matrix W of the eigenvectors of the matrix M gives the Fourier transform of the density 
modulation associated with the jth eigenvalue of M. We can thus produce an animation of the motion of the density 
modulation with the frequency tj and wave vector k by computing the time-dependent density 

6n (r, t) ~ e-^ l J2 e i(k+G) ' r 'Fjv,iv (k + G) W GJ (k) (23) 

G 

at several values of t in one period of the motion and then superimposing this density to the ground state density 
given by the Hartree-Fock calculation (eq. IjlipV We thus compute n (r. t) — Ubf (r) + aSn (r, t) where a is choosen 
sufficiently small to satisfy \Sn (r, t)\ < uhf (r). 



E. Orbital magnetisation and magnetic susceptibility 

The orbital magnetisation and susceptibility can be obtained from the dependence of the Hartree-Fock energy with 
filling factor. At T — OK, the orbital magnetisation per electron is given by 

1 fdE\ , s 



A V<9_B 

while the orbital magnetic susceptibility per particle is 

dm 



where E is the ground state energy of the 2DEG. 

When all Landau levels are taken into account, the Hartree-Fock energy per particle, excluding the Zeeman contri- 
bution, is (for < v* < 1) 



E T H °/ AL 
N T 



M<N a v 7 v 7 

4; (5) EE E w<G = o, 

v 7 a M<NM'<N 
v 7 M<N 

) £ [iW (G) (1 - fc, ) - Xn,n (G)] I (p (G)) | s 



(26) 



2^ Ui / 



G 



where 



min(7V,M) I 2 



Jo (qt x y) , (27) 



is the exchange interaction between electrons in Landau levels A and M, and a is the spin index. To derive this 
formula, we have assumed that all levels below A are occupied and that only the spin level a = — 1 in Landau level 
N is partially occupied. The total number of particles Nt is assumed fixed. 
With the definitions 

(28) 
(29) 
(30) 

(31) 

G 



A (A) 






M<N v 7 


Ai (A) 


= £ £ ^m,m> (g 




M<JV M'<iV 


A 2 (A) 


= £ XM ' N ( G = °) ' 




M<iV 


F (A, v*) 
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Eq. Il26[l can be written as 



E T H ° F TAL 



1 



v \ k£ 



1\ v* 
N + - I — ftw c 
v 

,2 



(32) 



- 3 ^ (iV) - - - A 2 (iV) + - - F (O , 



where 



F ("* ) = 2^ E ^ (G) (1 - J G , ) - Xjv.jv (G)] | (p (G)) | 

G 



(33) 



is the (dimensionless) Hartree-Fock energy per particle in the partially filled level i.e. the energy given in Eq. (|12J) . 

Differentiating Eq. (|32[1 once (twice) with respect to the magnetic field at constant density and adding the spin 
contribution, we get the total magnetization (total spin susceptibility). We write 



A* = MB ( ) (Ml + M2 + M3 + M4) > 

Vm* / 



(34) 



where [Ib = eh/2mc is the Bohr magneton. The four contributions to the magnetization are defined (for < v* < 1) 
by 



Mi 

M2 
M3 



-8A - - 2 ( AT 



2-111- 



1 

" 2 
4A\ 



f - 



47V 



(35) 



= -^=r s Ax (N) — v^r s A 2 (AT) 
v2 



--. + 3JV 



--„ + 3iV 



4=FK). 



Similarly, we write (for < v* < 1) the total magnetic susceptibility as the sum of the four contributions 



47mm* c 



(Xi +X2 + X3 + Xi) 



(36) 



where 



Xi 

X2 

X3 
Xa 



-8A (N) + 8N N + 



1 



(37) 



-2iV 5 * 



■m 



3r — 1 



j,3/2 



3A^ 



= V2r 



~ (*/ + 6AT) {y*) - (u + 2N) ^ 2 ^^ - (u - 2N) ^ 5/2 ^p 
4 ov* ov * z 



The four contributions to the magnetization and susceptibility come respectively from the kinetic energy (fix, xi), the 
Zeeman energy (/i2,X2)> the exchange energy between electrons in the filled levels or between electrons in the filled 
levels and those in level A" (fJ,3,Xs) an d the Hartree-Fock energy of electrons in level AT (n4,Xi)- I n Eq. (|37|l . we 
have introduced the gas parameter r s = I / ima B where as — h 2 n/m*e 2 is the effective Bohr radius and the effective 
g factor g* . We remark that our units for the susceptibility can also be written as e 2 /Aimm*c 2 = [i* B jB v= \ where 
H* B = eh/2m*c and B v= \ is the magnetic field required to get a filling factor of v = 1 for a total density of n. 



III. NUMERICAL RESULTS 



We now discuss our numerical results for the physical quantities introduced above. 
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FIG. 1: Hartree-Fock energy per particle as a function of the filling factor v* of the partially filled level for the bubble and 
stripe phases in higher Landau level. The inset in each figure indicates the critical filling factor Vi,j of the partially filled level 
at which the transition between bubble states or between the last bubble state and the stripe state occurs. 



A. Hartree-Fock energy 



Figures l(a)-(c) shows the Hartree-Fock energy per electron as a function of filling factor for the bubble and stripe 
phases in Landau levels N = 2, 3and 5. The filling factors at which the transitions between states occur are indicated 
in an inset in each figure. In all cases, the transitions are of first order and the transition scenario is that described 
in our introduction: the ground state evolves from the Wigner crystal at low filling, through a succession of bubble 
states with increasing values of guiding centers per bubble, and finally into the stripe state near half filling. The 
number of possible bubble states increases with Landau level index N. Our results for N = 2, 3are in agreement with 
previously obtained Hartree-Fock energies for the bubble states^. Energies for N — 5 have not been computed before. 
As shown in Ref. 0, the Hartree-Fock results are similar to those obtained by a DMRG calculation except that the 
region where the stripe phase is stable is wider in the DMRG calculation and, consequently, the bubble state with 
N + lelectron per bubble disappears. Oscillations in the Laguerre polynomial entering the form factor in Eq. © 
make it difficult to compute numerically accurate phase diagram for higher Landau level indices using our equation 
of motion technique. 



B. Density pattern 



Figure 2 shows how the density pattern given by Eq. (fill) evolves as the filling factor increases in Landau level 
N = 3. When the bubbles are relatively far apart, the density pattern on each crystal site is well approximated by 
the trial wave function of Fogler and Koulako\£. For a bubble of M electrons in Landau level N, that wave function 
is the Slater determinant 



*at (ri,r 2 , ...,r M ) 



^AM) 0i) 
<Pn,i (ri) 



Vjv.o (1*2) 
<Pn,i (1*2) 



<Pn,o (rjif) 
<Pn,i (rjf) 



<Pn,m-i ( r i) <Pn,m-i (ra) ••• <Pn,m-i ( t m) 



(38) 



,„L„. ,„ C_\ _ r» (r\\ m - N \ -r 2 /4£ 2 c i(N-m)e T \m-N\ 

wnere ipN,m [ r ) = ^N,m [j) e ' e L (N+ m -\m- 

tron in the symmetric gauge A 



(N+m-\m-N\)/2 

(^)is the normalised wave function of an elec- 
(— J3o2//2, Bqx/2,0) with Landau level ./V and angular momentum m. The one- 
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FIG. 2: Single-particle density of states and density pattern for bubble state with M = 1,2,3,4 electrons per bubble and for 
the stripe states in Landau level N = 3. The filling factors are (a) v* = 0.113; (b) v* — 0.173;(c) v* = 0.320; (d) v* = 0.406; 
(e) v* = 0.430. (See fig2.png file). 



particle density 

|4V(r,r 2 ,...,r M )| 2 , (39) 
is just 

m=M-l 

n N {r)= J2 WN, m {v)\ 2 - (40) 

m— 

As the filling factor increases in the Mth bubble state, the outer ring from each bubble gets closer to its neighbors and 
the density pattern given by Eq. (|40|l is strongly modified. When rings from two adjacent bubbles start to overlap, 
there is a transition to the M + 1th bubble state. In Landau level N, the last bubble state before the transition 
to the stripe phase has M = N + 1. We note that the stripe state obtained in our Hartree-Fock approximation is 
not the quantum Hall smectic state. It has density modulations along the stripes and can be described as an highly 
anisotropic Wigner crystal. However, it is likely that this state is unstable due to quantum fluctuations to the smectic 
stated. 



n N (r) = 



n 



dv, 



C. Density of states 

Figure 2 shows the behavior of the single-particle density of states g (E) defined in Eq. (|16fl for Landau level N = 3 
and for filling factors corresponding to bubble states with M = 1,2,3 and 4 electrons per bubble. In Landau level 
N = where the Wigner crystal state with M — lis the ground state for all filling factor, our calculation for g (E) 
reproduces the well-known Hofstadter butterfly structure. For v* = q/p, the density of states consists of p subbands of 
which q low-energy bands are filled and separated by a gap from the remaining p — q subbands. In Fig. 2, we see that 
the density of states in the bubble state has a different structure. The number of low-energy subbands that can be 
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resolved by our numerical calculation clearly corresponds to the number M of electrons per bubble. A measurement 
of the density of states, by a double- well tunneling experiment or photoluminescence, could serve to determine M. In 
computing the tunneling current, one would have to take into account other processes involved, particularly phonon 
shake-up. This is however beyond the scope of our paper and is left for future research. We note for now, however, 
that the typical phonon bandwidth (discussed below) is not larger than the splittings among the occupied bands in 
the density of states, so that the multiplicities of peaks associated with different bubble states are likely to survive 
such effects. 



D. Magnetic susceptibility 

The Hartree-Fock energy shown in Fig. 1 has discontinuous change of slope and curvature at the transitions 
between bubble states with different number M of electrons per bubble. This, in turn, gives rise to discontinuities in 
the behavior of the magnetic susceptibility with filling factor. To show this, we compute the susceptibility in Landau 
level N = 3 by numerically evaluating the first and second derivatives of the function F (y*) defined in Eq. Q33[) . We 
use the parameters Ao (3) = 9/2, Aj (3) = 6.29 and A2 (3) = 1.32 appropriate to N = 3 and assume a typical total 
electronic density of n = 3 x 10 n e/cm 2 so that the gas parameter is r s = 1.011. The other parameters are g* = 0.45 
and the effective mass m* = 0.067m. 

Figures 3(a) and 3(b) show the behavior of the four contributions to the magnetic moment and magnetic susceptibil- 
ity defined in Eq. (|35|l and Eq. I|37|). The singular behavior (contributions 114 and \i) comes, in both cases, from the 
Coulomb interaction between electrons in the partially filled level. The magnetic moment and susceptibility change 
discontinuously at the transition between bubble states and between bubble and stripe state. In both cases the effect 
is sizeable. Consequently, the transitions between bubble states should, in theory, be observable in a measurement of 
the magnetization or magnetic susceptibility^. 

E. Collective excitations 

We now discuss the collective excitations of the bubble states. We compute the density-density response function 
of Eq. i|22|) and follow its poles when the wave vector k is varied along the edges of the irreducible Brillouin zone 
of the triangular lattice of the crystal. Figures 4(a) to 4(f) show examples of the dispersion relation we obtain. In 
these figures, the dispersion is plotted along the path r — J — X — V corresponding to the wavevectors (k x ,k y ) = 
(0, 0), v(^?3' l)' ^"(773' 0) wri ere a is the lattice constant. The wavevector k represents the total distance, in 

reciprocal space and in units of 27r/a where a is the lattice constant, along the path r — J — X — V from the origin 
r. For a given value of k, Eq. I)32|l provides us with a way to image the motion of the density in a particular mode. 

We remark that our calculation does not include disorder so that the lowest-energy mode in figs. 4 is a gapless 
phonon mode. For small wave vector, this phonon mode has the typical uo ~ fc 3 / 2 dispersion of a Wigner crystal. In 
the presence of disorder, this phonon mode would be gapped. 

Upon entering a state with M electrons per bubble, the dispersion and maximal frequency of the phonon mode first 
increases with the filling fraction. Near the critical filling factors v* M ^ M+l computed above in the HFA, the phonon 
mode starts to soften at a finite wave vector k. This indicates that the bubble states become locally unstable when 
the outer rings of adjacent bubbles start to touch. For the cases = 2,3,5 that we have studied, we find that the 
phonon mode frequency vanishes at this wavevector for v* >v* M ^ M+1 so that a second order transition is preempted 
by a first order one. This is also the case for the transition between the M = N + 1 bubble state and the stripe state 
as we show in Fig. 5 where we plot the dispersion relation of phonon mode of the stripe state for wave vector k along 
the direction of the stripes (k x = n/t; where £ is the interstripe distance). The dispersion is shown for several values 
of the filling factor from v* = 0.42 to v* — 0.50 in Landau level N = 2. It is clear that the dispersion of the phonon 
(and in particular the region around the roton miminum) does not change significatively near v = 0.428, the filling 
factor at which the Hartree-Fock calculation predicts a transition between the M — 3 bubble state and the stripe 
state. 

For Landau level index N = 1, we find an interesting exception in that the softening of the phonon mode occurs 
for v* < v*_^2- This indicates the possible existence of some new charge ordered state. We did not study this case 
further, however, since our approximations do not include correlations responsible for the correlated liquid states that 
are present in N = 1—. 

Another indication of the limitations of the HFA and TDHFA is that the last bubble state (M = N+ 1) is stable, in 
contrast with the DMRG result. The DMRG method includes correlations neglected in the HFA and is, in principle, 
more exact than the HFA. In any case, for values of v where a bubble state is stable, we expect the TDHFA to give 
a reasonable approximation of the dispersion of its collective excitations. 
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FIG. 3: (a) Change in the magnetic moment and (b) magnetic susceptibility with partial filling factor v* in Landau level N=3. 
The four contributions i = 1 (solid), 2 (dashed), 3 (dashed dot), and 4 (squares) to /i or x are defined in eqs. H35H and 13711 . 

We see from figs. 4 that, above the phonon mode and in the region where the M = lbubble state is stable, 
there is a series of higher-energy modes that are almost dispersionless except near filling factors where the transitions 
between bubble states do occur. As we explained in Ref. [lj, more and more of these modes appear in the response 
function when we increase the size of the matrix Mg,g (k)in the equation of motion (Eq. (|18fl ) for the response 
function XG,G' (k,w). As more modes appear, the previous ones are not shifted in energy so that these modes are 
not numerical artifacts associated with the truncation of the infinite dimensional matrix Mg,g (k). For higher values 
of M, a general trend of the dispersion (except close to transition points) is that there are M low-energy dispersive 
modes close in energy and a number of higher-energy much less dispersive modes. For large values of M, each bubble 
has a complex pattern of density modulations and a large number of reciprocal lattice vectors are needed to describe 
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FIG. 4: Dispersion relations of the phonon and higher-energy modes in Landau level N = 3 along the edges (path V — J — X — T) 
of the irreducible Brillouin zone of the triangular lattice. The filling factors and number of guiding centers per bubble are (a) 
v* = 0.106, M = 1; (b) v* = 0.161, M = 1; (c) v* = 0.167, M = 2; (d) v* = 0.206, M = 2; (e) i/* = 0.264, M = 3; (f) 
v* = 0.380, M = 4. 



its structure adequately. It follows that our numerical procedure is less precise in this case as can be seen from Fig. 
4(f). 

The small dispersion of the higher-energy modes suggests that they can be identified as local oscillations of the 
density within a bubble. This is indeed what can be seen from an animation of these modes using Eq. (|32l) . Figs. 
6 show several snapshots of these modes representing the motion of the density in the second, third, and fourth 
modes for M — 1 in Landau level N = 2. From these figures, we see that the higher-energy modes are density waves 
propagating along the rings of the bubbles. These modes can be readily identified by the number of wavelengths of 
these waves enclosed in the perimeter of the rings. For example, the second mode corresponds to the case where 
the perimeter of the rings enclosed two wavelengths, the third mode to the case where the perimeter enclosed three 
wavelengths and so on. When the separation between the bubbles decreases in a given M bubble state, the density 
waves on each bubble become more and more coupled and the dispersion of the higher-energy modes can be very 
pronounced. 

The frequency of the higher-energy collective modes change discontinuously with the transition between the bubble 
states. A measurement of the frequency of these modes could then provide yet another signature of the bubble 
states. It would be interesting if variations in the gap w (k — > 0)of these higher-energy modes could be detected 
experimentally. Recently, Lewis et al^i reported the observation of an absorption peak in a measurement of the 
microwave absorption of a ultrahigh mobility two-dimensional electron system in Landau levels N = 2 and N — 3. A 
peak of absorption occurs at a frequency that decreases as the center (v* — 0.5) of the Landau level is approached 
from cither above or below. The resonances in the absorption are sharpest for filling factors v* around 1/4 and 3/4 
corresponding to the values where the RIQHE occur in transport experiments 2 . At these filling factors, the resonance 
frequency is approximately / = 500MHz. It is natural to associate this resonance with a pinned phonon mode of the 
bubble states (of either electrons or holes for v* < 0.5or v* > 0.5). 

The higher-energy modes of the bubble states are gapped at w (k = Q) = 0. In Fig. 7, we plot the gap / (k = 0) = 
ui (k = 0) /2ir in the first high-energy mode as a function of filling factor v* for Landau levels N = 2, 3, 5. To convert 
our frequencies (in units of e 2 /Hk£) to GHz, we assume k — 12.9 for the dielectric constant of the host semiconductor 
(GaAs) and take a typical density of n = 3.2 x 10 11 electrons/cm 2 so that e 2 jhnl — 3827/ y^GHz where v is the total 
filling factor of the 2DEG. At filling factor v, the magnetic field is given by B — nh/ev = li.2/v Tesla. As Fig. 7 
shows, the transitions between bubble states lead to abrupt changes in the gap frequency when the filling factor is 
varied. The discontinuities are more pronounced for lower Landau level indices N . The frequency range of the first 



12 



0.050 - 



0.000 







► 

► 


* \ 


► 

► 

► T 

• V" 


- * v 3* 


* 
* 


■ v*=0.42 
a v*=0.43 
▼ v*=0.44 

v*=0.50 











0.1 



0.2 



0.3 



0.4 



0.5 



k /(27i/a) 



FIG. 5: Dispersion relation of the phonon mode of the modulated stripe state in Landau level N — 2 for wave vector k along 
the direction of the stripes (k x = n/£). Here, a is the period of the modulations along the stripes and £ is the interstripe 
distance. 



higher-energy mode is in a range of frequencies that can be obtained in a microwave absorption experiment; thus, such 
experiments can in principle proble the transitions among bubble states. We remark that other effects not included 
in our calculation, particularly the modeling of a finite width for the electronic wave function in the quantum well, 
could lead to a reduction of the gap frequency. This would make behavior such as that seen in Fig. 7 easier to detect. 

To complete our calculation, it would be necessary to find out whether these higher-energy collective excitations 
are really detectable by a measurement of the absorption power. For this, one needs to compute the longitudinal 
conductivity a xx 'm the presence of the disorder potential since otherwise (by Kohn's theorem) only the cyclotron 
mode will show up in a xx . This calculation is beyond the scope of this paper but is in progress and will be reported 
elsewhere. 



IV. CONCLUSION 



In this work, we have studied several physical properties of the bubble states that form in higher Landau levels. 
We have computed the energy and density pattern of the 2DEG ground state, in Landau levels N = 2, 3, and 5, 
when the filling factor of the partially filled level is gradually increased from v* = to v* = 0.5. In the Hartree-Fock 
approximation, the Wigner crystal at small v* evolves into the modulated stripe state near v* = 0.5 by passing through 
a succession of bubble states with increasing number of guiding centers per bubble. In all cases (N > 1) that we have 
studied, the transitions are first order. 

We have shown that several physical quantities such as the single-particle density of states, the magnetization, the 
magnetic susceptibility and the dispersion of the collective excitations change discontinuously at the transition between 
these different ground states. We believe that these abrupt changes can be detected experimentally. In particular, we 
noticed that the density of states structure in the bubble states has features that allow one to determine, at least in 
principle, the number of guiding centers per bubble. 

We have studied more closely the collective excitations of the bubble states and, in particular, the structure of the 
higher-energy modes (i.e. those modes above the phonon mode). In light of the recent microwave experiments by 
Lewis et al»ii, we think that some of these modes are likely to be accessible experimentally. More work is needed to 
compute the real weight of these modes in the absorption. If these higher-energy modes can be detected in microwave 
absorption experiments, they will show a discontinuous change of the frequency / (k = 0) at each transition between 
the bubble states. 



13 




FIG. 6: Snapshots of the motion of the density in the second, third, and fourth higher-energy collective modes of the M = 1 
bubble state in Landau level N — 2. The filling factor is v* = 0.113. (See fig6.png file) 
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FIG. 7: Frequency / (k = 0) of the first higher-energy collective mode of the bubble states versus filling factor v* in Landau 
levels N = 2, 3, 5. We have assumed a density of n = 3.2 x 10 11 e/cm 2 . 
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